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Abstract. 

A multitransitional study with the BIMA interferometric array was carried out toward the starless core found 
in the L673 region, in order to study the small-size structure of the cores detected with previou s sing l e-dis h 
observations, which provides us with a test of the predictions of the chemical model of Taylor et al. (1996; 1998). 
We detected emission in the CS (J=2-^l), N2H^ {J=1~*0), and HCO^ (J=1^0) lines. Several clumps of size 
^ 0.08 pc wer e foun d for each line distributed all over the region where previous single-dish emission was found 
(Morata et al. 1997). Each molecular transition traces differently the clump distribution, although in some cases 
the detected clumps are coincident. The distribution of the N2H+ emission and the single-dish NH3 emission are 
coincident and compatible with an origin in the same gas. The large fraction of missing flux measured for the CS 
(2-^1) transition can be explained if the cloud is formed by a clumpy and heterogeneous medium. Four positions 
were selected to derive the abundance ratios [N2H^/CS] and [HCO'''/CS] from the molecular column density 
determinations, and to compare them with the values predicted by the chemical model. The model was able to 
explain the interferometric observations, and, in particular, the chemical differentiation of the detected clumps and 
the coincidence of the NH3 and N2H+ emissions. The lack of HCO"'' towards the two selected positions that trace 
the more evolved clumps cannot be accounted for by the model, but it is possibly due to strong self-absorption. 
We propose a classification of the studied clumps according to the stage of chemical evolution indicated by the 
molecular abundances. 
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1. Introduction 

Low-mass star formation takes place in dense cores of 
molecular clouds (Beichman et al. 1986; Benson and Myers 
|1989| ). The emission of several molecules, such as CS, 
NH3, and HCO"^, known to be good tracers of high den- 
sity molecular gas, has been used to study these dense 
cores. However, very soon large discrepancies between the 
emission of these molecules were found in some sources 



(Zhou et al. |1989| Myers et al. |l99l| ). To clarify the in- 
trinsic differences between the emission of these molecules 
and how this emission is related to the actual distribution 
of the high density gas, we began a systematic compari- 
son between the emission of the CS ( J=l— *0) and NH3 
{J, K)={1,1) transitions unde r simi lar conditions of an - 
gular resolution (Pastor et al. 1991 ; Morata et al. 1997). 



The comparison of the distribution of the CS and NH3 
emission in 14 condensations of 11 star- forming regions 
confirmed the discrepancies. In particular, we found that 
there is a separation ~ 0.2 pc between the emission peaks 
of both molecules; regions traced by CS are larger than 



those traced by NH3; and CS lines are ~ 0.5 km s ^ wider 
than those of NH3 . 

To explain these results, we developed a chemical 
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model (Taylor et al. 1996) in which high density conden- 
sations, or dense cores, are formed by clumps ^ 0.1 pc in 
size, which would be unresolved at moderate angular reso- 
lution observations such as the ones made in our study, of 
different mass, age, size and density. Most of the clumps 
disperse before NH3 abundances build up to significant 
levels, though these clumps contain substantial CS, so its 
emission should be observable. A few clumps, those suf- 
ficiently long lived, or being in a more advanced stage of 
physical and chemical evolution because of being denser 
or more massive, form a significant content of NH3, while 
CS abundance decreases with time. These clumps would 
possibly continue their evolution to eventually form stars. 
This model would account for the difference in size and 
separation between emission peaks of the CS and NH3 
molecules. 

Our chemical model explains the differences as a re- 
sult of the speed at which the molecules form. Therefore, 
a classification could be made between "early-time" 
molecules and "late-time" molecules, according to the 
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Table 1. Transitions observed 









RMS 


Spectral 




Molecule 


Transition 


V 


per channel 


resolution 


Detected? 






(GHz) 


(Jy/beam) 


( km s"^) 




HCO+ 


1-^0 


89.18852 


0.17 


0.33 


Yes 


C3S 


16^15 


92.48849 


0.18 


0.32 


No 


N2H+ 


1^0 


93.17402 


0.16 


0.32 


Yes 




07 — > ' 6 


yo.o / uiu 


n 1 Si 

U. lo 


U.Ol 


IVn 


ocs 


8^7 


97.30121 


0.22 


0.30 


No 


CS 


2^1 


97.98097 


0.18 


0.30 


Yes 


HNCS 


9l,8^8l.7 


105.74386 


0.25 


0.28 


No 


C2S 


98^87 


106.34774 


0.24 


0.28 


No 


HC3N 


12-»11 


109.17363 


0.26 


0.27 


No 


SO 


23^12 


109.25218 


0.25 


0.27 


marginally 


OCS 


9-^8 


109.46306 


0.25 


0.27 


No 



time at which these molecules reached their peak frac- 
tional abundance. We examined whether there were other 
potentially observable molecules that should show ex- 
tended emission like CS, or more compact emission as 



NH3 (Taylor et al. 1998 ). Several candidates were found 
for both groups: e.g. HCN, H2CO or HC3N are found in 
the CS family, whilst HCO+, SO, NO or N2H+ are in the 
NH3 family 

In order to study the small-size structure of the cores 
detected with the single-dish observations, which also pro- 
vides us with a test of the predictions of the chemical 
model, we carried out high angular resolution observa- 
tions towards L673 of several transitions corresponding 
to molecules of both families of species. The L673 high 
density condensation was previously detected in the CS 
(J=l— s-0) transition at moderate angular resolution by 
Morata et al. (1997). The L673 CS condensation is located 
at a(J2000) = 19'^20"5X1, (5(J2000) = -hll°15'29'.'5, at 
~ 9f5 to the southeast of IRAS 19180-1-1116. This conden- 
sation is well suited to our purposes for several reasons. A 
high density condensation has been detected nearby in th e 
NH3 (J,A') = (1,1) transition by Sepiilveda et al. (|2002|) , 
which shows the same general behavior found in the other 
sources of our survey: a separation from the CS emission 
peak ~ 2', and a NH3 condensation size clearly smaller 
than the CS condensation size. It is located nearby, the 
estimated distance to the L673 cloud is ^ 300 pc (Herbig 



and Jones, 1983), which helps in the study of the smaller 
size structure of the cloud. Finally, the condensations de- 
tected in CS and NH3 show no signs of tracers of star 
formation, such as infrared or radio continuum sources, 
Herbig-Haro objects or molecular outflows, which indi- 
cates that it is probably a quiescent core, maybe in the 
first stages of collapse, before forming a stellar core or a 
Class star, and thus fulfilling the conditions of the chem- 
ical model we developed. 



2. Observations 

The observations of the L673 region were carried out 
in 1998 May with the 10-antenna BIMA arrayQ at the 
Hat Creek Radio Observatory in the C configuration. 
The phase calibrators were 1751+096 and 3C395. In or- 
der to include the positions of the CS (J=1^0) and 
NH3 (J, Ar)=(l,l) transitions singlc-dish emission peaks 
of t he m aps by Morata et al. (1997) and Sepiilveda et 
al. ( 2002 ) , we made a two-point mosaic with each of the 
fields centered approximately at the position of these two 
peaks. Thus, one of the fields was centered at the posi- 
tion a(2000) = 19'>20'°52!2, and (5(2000) = +11°13'57", 
and the other located 90" to the North. Three frequency 
setups were used, centered at 91, 96 and 107 GHz. The 
digital correlator was configured to observe simultaneously 
several molecular line transitions at moderately high spec- 
tral resolution, ^ 0.3 km s^^. The target molecular lines 
were HCO+ (J=1^0), N2H+ (J=1^0), CS (J=2^1), 
SO (Jk = 23 I2) and C^^O (J=l->0). System tem- 
peratures for the 91, 96 and 107 GHz setups were in the 
180-500 K, 250-700 K and 500-1000 K range, respectively. 
The calibration and data reduction were performed using 
the MIRIAD software package (Sault, Teuben & Wright 
1995). Mosaic maps were done with the visibility data 



weighted by the associated system temperatures and using 
natural weighting, and applying the primary beam correc- 
tion. The resulting synthesized beam for the 91, 96 and 
107 GHz setups were 12'.'7 x 9'.'6, PA = -14°; 13'.'3 x 10'.'5, 
PA = -17°; and IS'.'l x 9'.'3, PA = -8°, respectively. The 
transitions of all the observed lines, their velocity resolu- 
tion and the achieved rms noise with this velocity reso- 
lution are listed in Table |l} For the continuum emission 
we used the data from the 96 GHz correlator setup, which 



^ The BIMA array is operated by the Berkeley-IUinois- 
Maryland Association with support from the National Science 
Foundation. 
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Fig. 1. Integrated emission of the CS ( J=2^1) line for the 
Vlsr range 5.6-9.8 km s'^. Contours are 0.96, 1.20, 1.44, 
1.68, 1.92, 2.16, 2.40, 2.64 Jy beam"! km s'^ The beam 
(20" X 15'.'8) is shown in the bottom left-hand corner. The 
upright crosses indicate the nominal position of the single- 
dish emission peak of the CS ( J=1^0) line (north) and 
NH3 {J,K) = {1,1) line (south). The tilted crosses indicate 
the position of the four points selected for further study 
(see Table |^). The bordering contour indicates the two 
fields observed 



provided a bandwidth of 800 MHz centered at the fre- 
quency of 95.9 GHz. No emission was detected at above 
8 mJy beam~^ (at a 3-cr level). 



3. Results 

3.1. Morphology of the molecular emission 

We detected emission in three of the transitions: CS 
(J=2^1), N2H+ ( J=1^0), and HCO+ (J=1^0). We de- 
tected marginally, at ^ 2(7 level, the SO (23^12) transi- 
tion. The remaining transitions of Table 0, which where 
observed mainly because they were found inside the fre- 
quency range of the observations, were not detected. It 
must be noted that all these undetected transitions have 
an upper level energy with temperatures higher than 15 K, 
and most of them, OCS is the exception, have also a high 
dipole moment. Thus, the temperature and density condi- 
tions of the L673 core, and the low abundances expected 
for these molecules make them very difficult to detect. 



Fig. 2. Same as in Fig. for the N2H+ (J=1^0) line 
for the Vlsr range —1.4-14.2 km s^^. Contours are 
2.55, 3.06, 3.57, 4.08, 4.59, 5.10, 5.61, 6.12, 6.63, 7.14 Jy 
beam"^ km s^^. The beam (20" x 14"6) is shown in the 
bottom left-hand corner. 



Figures |l| to || show the maps of the zero-order mo- 
ment (integrated intensity) of the emission of the detected 
molecules. These maps were obtained after convolving the 
original maps with a Gaussian function, with a resulting 
beam of ^ 20" x 15", in order to make a more meaningful 
comparison between the emission of each molecule. 

A clumpy distribution of the emission is shown in the 
integrated intensity map of the CS (J=2— >1) transition 
(Fig. ^) . The more intense emission is found in two clumps 
near the position of the single-dish NH3 {J,K)—{1,1) 
emission peak: one almost coinciding with its nominal po- 
sition, and the other peaking ^ 30" to the west. Weaker 
emission is found in a clump near the position of the single- 
dish CS (J=1^0) emission maximum, ~ to the north. 

Figure]^ shows the integrated intensity map of N2H+ 
(J=l— >0). Two clumps of emission are found, one al- 
most coincident with the nominal position of the NH3 
{J, K) = (l,l) emission peak, and the other, more intense, 
~ O'.Q to the south. No emission is detected around the 
position of the CS ( J=1^0) single dish emission peak. 

The integrated intensity map of the HCO+ (J=l— >0) 
emission (Fig. |^) shows several isolated emission enhance- 
ments distributed around the positions of the CS and 
NH3 single-dish emission peaks. The most intense emis- 
sion tends to be around the position of the CS (J=l— >0) 
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Fig. 3. Same as in Fig. 0for the HCO+ (J=1^0) line for 
the Vlsr range 5.6-9.8 km s~^. Contours are 0.96, 1.08, 
1.20, 1.32, 1.44, 1.56, 1.68 Jy beam"! km s'^ . The beam 
(20" X IS'.'l) is shown in the bottom left-hand corner. 



emission peak or in-between the two peaks, with a N-S 
elongation. Another intense enhancement is found ^ I'.l 
south-east of the NH3 emission peak nominal position. 

Comparing the integrated intensity maps of the three 
molecules, we found that N2H"'" emission was more con- 
centrated than CS and HCO"*" emission, which were found 
more spread all over the region. CS (J=2— >1) emission in 
the southern region coincided with the northern clump 
of N2H"'' emission. Moreover, the northern N2H+ emis- 
sion peak was found between the two CS ( J=2-^l) local 
emission enhancements. CS (J=2^1) and HCO+ emis- 
sion were found to coincide closely at the northern part of 
the region, especially around the CS (J=l— *0) emission 
peak. On the contrary, N2H+ and HCO+ emission coin- 
cided only marginally in general. Finally, we were also 
able to find an emission enhancement in the marginally 
detected SO emission very near the position of the N2H+ 
northern clump. 

Table ^ lists the four positions we selected in the 
mapped region in order to compare the physical parame- 
ters of the gas traced by each molecule. The selected po- 
sitions were related to the most prominent higher resolu- 
tion clumps found in our observations, and corresponded 
to local intensity peaks in the emission of the detected 
molecules. The four positions were labeled as S, E, W 




5 10 15 5 10 15 5 10 15 5 10 15 
Velocity (km/s) 

Fig. 4. Spectra of the CS (J=2^1), HCO+ (J=1^0), 
and N2H+ ( J=1^0) transitions at the four selected posi- 
tions in our mapped region. 



and N, and roughly correspond to the South peak posi- 
tion of the N2H+ integrated intensity map, the East CS 
local maximum, the West CS maximum, and the North en- 
hancement in the HCO+ emission near the CS ( J=l— >0) 
emission peak, respectively. 

3.2. Kinematic structure 

Figure ^ shows the spectra obtained for the three tran- 
sitions detected in the region in the four selected points. 
Table || lists the fitted line parameters for each detected 
transition. CS emission is detected, with varying inten- 
sity, in all four positions, whereas HCO+ emission is only 
detected at the N and W positions, and N2H+ emission 
at the S and E positions. Line center velocities for CS 
and HCO+ at the positions where a reliable fit could 
be obtained are compatible with being originated in the 
same bulk of gas. The difference in the line center ve- 
locity is ~ 0.3-0.5 km s~^, which is the velocity reso- 
lution of our observations. However, line center veloci- 
ties for the N2H"'' spectra were shifted > 0.7 km s~^ 
from the CS line velocities. Using the much better de- 
termined frequency for the main N2H+ (,/=1^0) compo- 
nent {Fi,F = 2,3 ^ 1,2) given by Caselh et al. ( |l995[) : 
93.17378 GHz, we found a good agreement with CS spec- 
tra line velocities, ^ 0.3 km s^^. We have corrected the 
N2H"'' line velocities accordingly, as Table ^ shows. Line- 
widths are < 0.8 km s~^, very similar to the observed 
values in the single-dish observations, except for the W 
clump where velocity dispersion is higher, 1.0 km s^^. 



3.3. Physical parameters 

We were able to estimate the excitation temperature for 
the N2H+ data at the two positions, S and E, where we 
had good signal-to-noise spectra. The values obtained at 
both positions are ~ 4 K (assuming a beam filling factor 
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Table 2. Selected positions 



Position map counterpart 

a (J2000) 5 (J2000) 



South (S) 


19''20™51.61' 


11°13'40'.'5 


N2H+ 


(J= 


1^0) south peak 


East (E) 


19'^20™52.29' 


11°14'10'.'5 


CS ( J-- 


=2^1) 


eastern peak 


West (W) 


19*^20'" 50.38' 


11°14'14'.'5 


CS {J-- 


=2- 


>1) 


western peak 


North (N) 


19''20™ 53.24' 


11°15'30'.'5 


HCO+ 


{J-- 


=1- 


-^0) NE enhancement 



Table 3. CS (J=2^1), HCO+ (J=l->0), and N2H+ (J=1^0) lines and physical parameters obtained at the four 
selected positions in our region 



Molecule & 


Pes. 


Tmb 




Vlsr 


AV 


T (*) 






Transition 




(K) 


(K) 


(kms"') 


(kms-') 






(10" cm-2) 


CS 


S 


0.35 ±0.12 




±7.4 ±0.1 


0.7 ± 0.2 




0.05-0.4 ("') 


12.0-26.9 


(7=2^1) 


E 


0.96 ±0.11 




±6.9±0.1 


0.5 ±0.1 




0.15-2.1 ('') 


25.3-84.8 




W 


0.73 ±0.10 




±6.8 ±0.1 


1.0 ±0.1 




0.11-1.1 ('') 


37.2-99.6 




N 


0.72 ±0.11 




±7.3±0.1 


0.8 ±0.1 




0.11-1.1 w 


27.1-72.2 


HCO+ 


S 


< O.37W 










< 0.39 


<4.73 


(7=1^0) 


E 


< 0.29(''' 










< 0.29 


< 3.60 




W 


0.48 ±0.13 




±7.2 ±0.1 


0.8 ±0.3 




0.07-0.6 ('') 


4.7-7.7 




N 


1.02 ±0.12 




±7.0 ±0.1 


0.6 ±0.1 




0.16-2.3 W 


7.7-19.1 


N2H+ 


S 




0.32 ± 0.05 


±7.0 ±0.1 


0.5 ±0.1 


0.25 ±0.10 


1.0 («' 


2.2-3.1 


(7=1^0) 


E 




0.24 ± 0.04 


±7.0 ±0.1 


0.5 ±0.1 (■'■' 


0.27 ±0.10 


1.0 («) 


1.7-2.4 




W 




< 0.20*''' 








< 0.7 W 


<2.0 




N 




<0.16(''' 








< 0.5 I'" 


< 1.6 



Derived from the transfer equation, where A = f[J{T^y^) — y(7bg)] is the "amphtude" (Pauls et al. 1983), / is the filling factor, the excitation 
temperature of the transition, the background radiation temperature, and J[T) the intensity in units of temperature. 

Optical depth of the main line. 
^''^ Beam averaged column density calculated from 

"* ^ 4^l^v^,'a expl-hv(J-l)/2kn,]-expl-hv(J+l)/2kT„.] J(T,,)-J(T^,) J 

where / is the beam filling factor, the dipole moment of the molecule, J the lower level of the transition, and J(T) the intensity in units of 
temperature. 

Derived from the transfer equation, assuming for the excitation temperature upper and lower values of 10 K and 4 K, respectively. 
Upper limit taken as 3a, where a is the sensitivity per channel. 
Intrinsic line width 
Total Une opacity 



/ = 1), which is compatible with the value obtained from 
the CS (J=1^0) observations, 4.2 K (Morata et al. |1997D , 
and slightly lower than the value obtained from th e NH 3 
{J,K)={1,1) observations, 5.7 K (Sepiilveda et al. |2002|) . 
The total line opacity of the N2H+ {J— 1^0) transition at 
both positions is ^ 1.0. We estimated an upper limit for 
the beam averaged column density using the value for the 
excitation temperature obtained from the CS (J=1^0) 
and N2H+ ( J=1^0) observations, Tex = 4 K, and assum- 



ing a filling- factor / = 1. However, the filling- factor could 
be / < 1. In this case, the excitation temperature would 
be greater than 4 K, but likely less than 10 K, because 
the kinetic temperature in low-mass dense cores is typi- 
cally about 10 K, and the excitation temperature will not 
be greater than this value. Thus, we used an upper value 
for the excitation temperature Tex = 10 K, which implies 
a lower limit for the beam average column density. The 
values for the line opacity, r, and the beam averaged col- 
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Fig. 5. Spectra of the CS (J^l^O) transitions obtained 
at the positions of the CS (J=1^0) and NH3 (J, K) = {1,1) 
single-dish emission peaks vifith the Yebes telescope (top), 
compared with the spectra of the CS (J=2^1) transition 
obtained at the same positions after convolving the BIMA 
data with a beam of 1.'9 equal to the single-dish beam 
(bottom). Note the different vertical scales. 



umn density, A^mb for the four selected positions are shown 
in Table |. 

3.4. Comparison with single-dish observations 

We compared the CS (.7=2^1) interferometric and 
(J=1^0) single-dish emissions. Figure || shows the single- 
dish CS (J=1^0) spectra obtained at the positions of the 
CS ( J=1^0) and NH3 (J, K)—{1,1) emission peaks, com- 
pared with the BIMA CS (J=2^1) spectra obtained at 
the same positions after convolving the data with a gaus- 
sian function to obtain a 1.'9 resulting beam, equal to the 
beamsize of the single-dish observations. Line center ve- 
locities agree within 0.3 km s~^, our spectral resolution in 
both sets of observations. Line widths are almost coinci- 
dent at the southern position, while the (J— 1—^0) line is 
~ 0.35 km s^^ wider at the northern position. In order to 
estimate the flux loss in our interferometric observations, 
we calculated the line intensity we should obtain for the 
(J=2^1) line from the (J=1^0) line opacity and exci- 
tation temperature obtained previously in the single-dish 
observations (Morata et al. 1997). We found that the mea- 
sured line intensity is ~ 8% in the northern position, and 
~ 12% in the south position from the estimated values. 
These results are compatible with the characteristics of 
the BIMA array, which is insensitive to structures larger 
than about 10 times the synthesized beamwidth (Wright 
1996|), which in our case is ~ l.'5-2' (0.13-0.17 pc), of the 
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Fig. 6. Contour map of the single-dish integrated line 
intensity of the CS (J=1^0) line (solid line) obtained 
with the Yebes telescope (Morata et al. 1997 ), and the 
single-dish main beam brightness temperature of the NII3 
(J, if) = (1,1) transition (dashed line) obtai ned with the 
Haystack telescope (Sepiilveda et al. 2002). The lowest 
CS contour level is 0.48 K km s^^, and the increment is 
0.05 K km s-i. The lowest NH3 contour level is 1.35 K, 
and the increment is 0.3 K. The crosses indicate the nom- 
inal position of the single-dish emission peaks of the CS 
(J=l->0) line (north) and NH3 {J,K)={1,1) line (south). 
The filled squares indicate the position of the four points 
selected for further study in the interferometric observa- 
tions 



However, the existence of such an extended compo- 
nent, responsible for ~ 90% of the single-dish line inten- 
sity, would not affect significantly the beam averaged col- 
umn densities determined from our interferometric obser- 
vations. Assuming a size of 120" for the extended com- 
ponent, its contribution to the beam averaged column 
density for a beam size of 17" would be of the order of 
90% X (17"/120")2 = 2%. Thus, the beam averaged col- 
umn densities determined from our observations are not 
significantly affected by the missing emission by the BIMA 
interferometer, assuming that it is due to an extended 
cloud component. 

Figure || shows the distribution of the single-dish CS 



(J=l-^0) (Morata et a l 
(Sepiilveda et al. |2002D 



1997) and 



NH3 
the 



order of the single-dish beamsize. 



(J,i^)=(l,l) 
emission m the same region 
mapped by the interferometric observations. We found 
that most of the strongest CS ( J=2^1) emission and par- 
ticularly all the detected N2H+ emission were enclosed 
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inside the half-maximum contours of the NH3 map, lo- 
cated at the southern region of the map. All four relative 
emission enhancements in these two transitions are also 
enclosed inside the beam of the single-dish observations. 
HCO"*" emission was fainter around this southern region, 
and marginal in the eastern, western and southern margins 
of the NH3 intensity contours. Around the northern re- 
gion, no N2H+ emission was detected, but CS and HCO+ 
emission were enclosed by the beam of the CS single-dish 
observations. Emission around this position was not as 
concentrated as around the southern region, but we found 
intense HCO+ and CS (J=2^1) emission very near the 
beam center. 



4. Discussion 

4.1. Structure of the core 

The integrated intensity maps of the CS (J=2^1), 
N2H+ (J=1^0), and HCO+ (J==1^0) transitions ob- 
tained with the BIMA interferometer showed that a much 
clumpier medium is revealed by the high angular reso- 
lution observations than b y the single-dish maps of CS 
(J=1^0) (Morat a et al. |l997|) and NH3 (J,if) = (l,l) 
(Sepiilveda et al. 2002 ) obtained with an angular reso- 
lution of ^ Several clumps of size ^ 0.08 pc are 
found distributed all over the region where the strongest 
emission of the single-dish observations was located, which 
would support the idea that at least part of the single-dish 
emission was originated in clumps of smaller size, < 0.1 pc 
(Taylor et al. 1996 ). Moreover, these small size clumps, al- 
though coinciding in some cases, are traced differently by 
each molecule, which would support the idea of a chemi- 
cal differentiation between each clump of gas depending on 
the local density and age. However, we found that ~ 90% 
of the single-dish CS emission is not detected by the inter- 
ferometric observations. As we have already pointed out, 
the BIMA array is insensitive to structures larger than 
~ 0.17 pc, which could mean that part of the emission 
could come from chemically young clumps of that size or 
larger, undetectable then by BIMA, but whose emission 
was detected by the single-dish observations. Another pos- 
sibility is that there could also be clumps smaller than the 
ones we found, of low density, and with emission not strong 
enough to be detected with the interferometer. 

In order to qualitatively test which of the aforemen- 
tioned two possibilities is more consistent with our ob- 
servations, we modeled the filtering effect of the BIMA 
interferometer response to a molecular core with an ho- 
mogeneous and extended geometry and a clumpy hetero- 
geneous geometry. The process of generating synthetic vis- 
ibility data matching that expected from true observations 
with the BIMA interferometer was performed by using 
the MIRIAD task UVGEN. This process is described by 



Girart & Acord (2001). Our synthetic observations had 
the same phase center as the L673 observations and were 
done for the BIMA C-array configuration. The total flux 
adopted for both cases, the extended and clumpy struc- 



tures, was similar to that expected for the CS (J=2^1) 
transition, calculated from the measured CS ( J^l^O) in- 
tensity. The results obtained were: 

1. Extended and homogeneous medium. The presence of 
the extended emission was modeled as single large 
component. We tested the filtering effect of the inter- 
ferometer for different sizes and different center posi- 
tions of the cloud with respect to the phase center. We 
found that in order to achieve a missing flux of about 
90%, the value estimated for the BIMA observations, 
large sizes were required: diameter ^ 2'. The resulting 
synthetic maps showed strong negative lobes extend- 
ing ~ 40" along the N-S direction, at both sides of 
the center of the synthetic cloud. This was due to a 
combination of the effects of source structure, with a 
large fraction of undetected flux by the interferome- 
ter, and visibility coverage of the BIMA interferome- 
ter for a source of low declination {6 — 11° for L673). 
These negative structures could not be removed in the 
cleaning process. The cleaned maps had absolute peak 
intensities of the negative structures with similar, or 
even higher, values than the positive structures. 
However the real BIMA CS (J=2^1) channel maps 
does not show such strong negative structures. At the 
central velocities of the CS line there are some negative 
lobes, but their intensity is only a ^ 40% of the positive 
peaks. 

2. Clumpy and heterogene ous medium . The chemical 
models of Taylor et al. ( 1996| ; 1998 ) assume that the 
clouds are formed by clumps of ^ 0.1 pc with different 
sizes, ages and masses. From this assumption it is ex- 
pected a cloud formed with only few massive clumps 
and several clumps less massive. So, we adopted a 
simplistic model of a cloud formed with one massive, 
thereby strong, clump and between 15 and 20 less mas- 
sive, thereby weaker, clumps. We found that even most 
of the emitting flux came from the smaller clumps, 
the BIMA synthetic resulting maps were clearly dom- 
inated by the strongest clump, and a large fraction of 
the emission from the weaker clumps was undetected 
by the interferometer. In some of the cases modeled, 
the total missing flux was of ^ 80% of the total flux. 
The residual negative structures of the synthetic maps 
were 30-50% of the positive. 

In summary, without going to an exhaustive analysis of 
the response of the BIMA interferometer to the different 
cloud structures, we find that in L673 the missing flux is 
most probably due to a clumpy heterogeneous structure. 

4.2. Chemical abundances 

In our maps, we found that all the detected N2H+ emis- 
sion was enclosed inside the half-maximum contours of the 
NH3 map of Sepiilveda et al. ( p002 ), with the strongest 
N2H"'" emission found inside the highest NH3 contours. 
The nominal position of the NH3 emission peak is placed 
between the two local N2H+ emission peaks, although the 
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Table 4. Column density ratios 



CQ 



Position 


[HCO+/CS] 


[N2H+/CS] 


N(CS)/N„,ax(CS) 


S 


< 0.21 


0.15 


30% 


E 


< 0.07 


0.05 


77% 


W 


0.10 


< 0.02 


100% 


N 


0.27 


< 0.02 


73% 



beam enclosed both maxima. The chemical model devel- 
oped by Taylor et al. ( |1996| ; |1998| ) marks the NH3 and 
N2H+ molecules as forming at late times, when sufficiently 
high densities are reached or in long lived clumps. Thus, 
they should be found in the same region. The coincidence 
between the distribution of the emission of N2H+ observed 
by the BIMA interferometer and the NH3 emission de- 
tected by single-dish observations, and the overall morpho- 
logical appearance of the NH3 and N2H+ emission seem to 
support these predictions. Moreover, both interferometric 
and single-dish emissions are highly concentrated, and are 
compatible with being originated in the same gas since 
the emission of N2H"'" could be coming from two clumps 
of size ~ 0.'6 X 1' enclosed inside the higher levels of the 
lower resolution emission of NH3. The SO molecule was 
also marked in our model as a late-time molecule. The 
near coincidence of the clump found in the marginally de- 
tected SO emission with the northern clump of the N2H+ 
emission would support that classification. 

Table |^ shows the column density ratios of N2H+ and 
HCO+ over CS for the four selected positions calculated 
from the column density determinations of Table ^, and 
the CS column density relative to the maximum value, 
found in position W. We observe that N2H+ abundance 
is ~ 5-15 % of that of CS for the S and E positions, and 
becomes lower at the N and W positions. [HCO^/CS] is 
^ 20% for the S, E, and W positions, and becomes ~ 30% 
at the N position, while there is a factor 3 variation in the 
CS column density. We can then differentiate between the 
N and W positions, where HCO"*" is clearly detected while 
N2H+ is not, and the E and S positions where N2H+ is 
clearly detected whereas HCO"*" is not detected. 



4.3. Modeling the chemical evolution 

We used the results of the model explained in Taylor et 
al. (1998) to compare the predicted abundances of HCO+ 
and N2H+ with respect to CS, with the values we ob- 
tained for the column density ratios, and see if some chem- 
ical age determination could be made for these positions. 
Figure ^ shows the predicted fractional abundances for 
several molecules of interest for values of collapse, B — 1, 
freeze-out, Fi? = 0.01, and final density n/ — SxlO** cm^'^ 
(see Taylor et al. 1996; 1998). We also show the molecular 
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Fig. 7. Left panel: Chemical fractional abundances (rela- 
tive to H nuclei) as a function of time for a free-fall col- 
lapse model halted at density nn = 5 x lO'' cm~'^. Initially 
riH = 1 X 10'^ cm~'^ and Ay — 0.5. Freeze-out parameter 
FR = 0.01, corresponding to an average value of the prod- 
uct of the dust to gas number density ratio and square of 
the grain radius (n^a^) = 2 x 10~^° cm~^. Right panel: 
Fractional abundances of selected species (relative to CS) 
as a function of time. 

abundance relative to CS for the HCO+, N2H+, and NH3 
molecules. 

We found that the abundance [HCO+/CS] can be ex- 
plained by the model in relatively short times, between 2- 
3x10^ years. Thus, a higher HCO+ fractional abundance 
would indicate a more chemically evolved clump. However, 
this does not seem to agree with the abundance of N2H''". 
First, N2H+ reaches high abundances at later times, and 
one would expect to find N2H"'" at more evolved clumps. 
Second, the model can not explain fractional abundances 
of N2H+ similar to those of HCO+ , and as high as those we 
found in the S and E positions. A slightly higher freeze-out 
parameter would provide higher relative N2H+ abundance 
with respect to CS, while not changing appreciably the 
fractional abundance, because CS is frozen-out to grains 
faster. In this case, we also obtain higher IICO+ relative 
abundances with CS. 

Figure ^ shows the results of the model modifying the 
final density value to 10 times the previous value, Uf = 
5 x 10^ cm^'^, and maintaining unchanged the values of 
B and FR. We find some differences with the previous 
results. CS maximum fractional abundance is similar to 
the previous model, while NH3, IICO+, and N2H''' reach 
slightly lower abundances. In any case, CS and NH3 show 
the same behavior found in the model of Taylor et al. 
( 19961 ). However, the higher final density produces a faster 
depletion of molecules, especially CS which has a steeper 
decline with time. 

At times ^ 2.5 x 10^ years, when CS has a high frac- 
tional abundance, and HCO+ has barely reached its max- 
imum fractional abundance, we find a relative abundance 
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Fig. 8. Left panel: Chemical fractional abundances (rela- 
tive to H nuclei) as a function of time for a free-fall col- 
lapse model halted at density tt-h = 5 x 10^ cm~^. Initially 
riH = 1 X 10'^ cm~^ and Ay = 0.5. Freeze-out parameter 
FR ~ 0.01, corresponding to an average value of the prod- 
uct of the dust to gas number density ratio and square of 
the grain radius {udO^) = 2 x 10~^° cm~^. Right panel: 
Fractional abundances of selected species (relative to CS) 
as a function of time. 



[HCO+/CS]~ 0.1, while N2H+ has still a low abundance. 
This could correspond to the W position of our map. At a 
later time, 3 x 10^ years, when N2H+ is still at low abun- 
dances, HCO+ abundance is practically the same while 
CS abundance is slightly lower, with a resulting higher 
[HCO+/CS] ratio. This could correspond to the N po- 
sition. At the time N2H+ reaches its peak abundance, 
~ 4 X 10^ years, HCO+ has diminished in abundance al- 
though it can still be observed. This would point to the S 
and E positions being at later stages of chemical evolution. 
Distinguishing between the two is difficult. The S position 
could be more evolved since it shows a lower CS column 
density, and a higher [N2H"'"/CS] ratio, which could corre- 
spond to the stage where N2H"*' is at higher abundances. 

However, even the modified model is not able to ex- 
plain the observed HCO+ low abundances obtained at po- 
sitions S and E with respect to the predicted values. The 
HCO+ molecule has a high dipole moment, and it has been 
found that the (J=l— >0) transition line profiles are of- 
ten self-absorbed in typical molecular cloud cores by fore- 
ground gas at lower density and temperature. The appar- 
ently sudden lack of HCO+ emission in those places where 
N2H+ or CS emission is intense may suggest that the 
HCO+ ( J=1^0) emission arising from the core might be 
absorbed efficiently by a cold low-density envelope around 
the core or a foreground cloud, as in the case of the NGC 
2264G region (Girart et al. po5o| ). 



4.4. Comparison with other studies 

Studies of chemical evolution have been made up to now 
with single-dish observations using several transitions of 
selected early- and late-time molecules, usually CS, N2H+, 
NH3, HCO+, SO, SO2, and HC3N. These studies of star- 
less cores show maps where emission seems to originate in 
several clumps distributed all over the region (Dickens et 
al. ^OOOD or in round-like clumps (Tafalla et al. ^0021 ). All 
these studies show evidences of chemical differentiation. 
Dickens et al. find that early- and late-time molecules are 
found at different positions and show differing behavior, 
i.e. early-time molecules are most abundant where late- 
time molecules are not, and viceversa. The authors explain 
their results invoking a difference in the chemical ages of 
different parts of their region. Tafalla et al. also find that 
their maps show a great coincidence between the emission 
of the N2H+ and NH3 molecules, while differing from the 
CS emission. 

Some of these studies also find evidence for differing 
depletion rates of molecules such as CS or NH3. Tafalla et 
al. find a central abundance drop of CS, a constant abun- 
dance of N2H"'', and a central enhancement of NH3 in 



their regions, while Bergin et al. (2001), in observations of 
the IC 5146 cloud, also detect molecular depletion of CS, 
while N2H+ remains in the gas phase with growing extinc- 
tion. Other numerical calculations, such as Aikawa et al. 



( 2001 ) also predict a depletion of certain species from the 
central regions. In our model, we would expect to find the 
more massive and chemically older clumps in the central 
regions. These clumps would have substantial quantities of 
NH3 and N2H+, while CS would have been depleted from 
the gas phase, which would show as a central depiction of 
CS. 

Thus, these results suggest that although our interfer- 
ometric observations lose information on the total emis- 
sion of the region, they are in agreement with the results 
obtained with single-dish telescopes, suggesting that in- 
terferometric observations are useful to study with high 
angular resolution the chemical evolution of inner struc- 
tures of starless cores. 



5. Summary 

We made a multitransitional study with the BIMA in- 
terferometric array of the starless core found in the L673 



(Morata et al. 1997), in order to test the chemical model 
of Taylor et al. (1996; 1998). The main results were: 



We detected emission in the CS(,/=2^1), N2II+ 
(J=1^0), and HCO+ (J=1^0) lines. We marginally 
detected emission in the SO {Jk = 2^ ^ I2) line. 
The high angular resolution interferometric observa- 
tions revealed a much clumpier medium than the lower 
resolution single-dish observations. Several clumps of 
size ^ 0.08 pc were found for each line distributed 
all over the region where the single-dish emission was 
found. 
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Each molecular transition traced differently the clump 
distribution, although in some cases the detected 
clumps were coincident. We found that the distribution 
of the N2H+ emission was completely enclosed by the 
half-maximum contours of the NH3 maps of Sepiilveda 
et al. ( ^002 ), and that both emissions were highly com- 
pact and compatible with being originated in the same 
gas. We also found a marginal coincidence with the SO 
molecule, which is also a late-time molecule. 
The BIMA interferometer detected only 9-12% of the 
CS (2— >1) emission. Modelling the filtering effect of 
BIMA, we found that a clumpy and heterogeneous 
medium could explain this effect. 
We estimated the abundance ratios [N2H+/CS] and 
[HCO^/CS] from the molecular column density deter- 
minations in four selected positions, labeled as S, E, W, 
and N, in order to compare them with the predicted 
values of the chemical model. We found that at the N 
and W positions there was a high HCO+ abundance 
relative to CS and a low N2H+ abundance, while at 
the S and E positions the case was reversed. 
The chemical model was able to explain the abundance 
ratios of [HCO+/CS] at the positions where HCO+ 
was detected, N and W. The best fit was for a model 
with a density at which collapse is halted oi nj — 5 x 
10^ cm^'^. In this case, we found that the predicted 
[N2H+/CS] abundance ratios were more in agreement 
with the observations. A high IICO+ abundance, but a 
low N2H+ abundance would represent an earlier stage, 
where the W and N positions could be found. The S 
and E positions would be in a later chemical stage on 
account of the higher N2H+ abundance, although a 
more precise differentiation is hard to make. 



Thus, the model described in Taylor et al. ( |l996t |1998D 
applied to a starless core with no apparent signs of other 
star formation tracers, such as molecular outflows or ra- 
dio continuum sources, was able to explain the interfer- 
ometric observations, if we take into account the likely 
self-absorption of HCO+. In particular, the chemical dif- 
ferentiation of the detected clumps and the coincidence 
of the NH3 and N2II+ emissions agree with the model. 
It also enabled us to explain the HCO+ and N2H+ emis- 
sions and to propose a classification of the studied clumps 
according to the stage of chemical evolution indicated by 
the molecular abundances. 

Further multitransitional observations with higher sen- 
sitivity would be necessary in order to determine with 
more detail the physical structure of the region and to pro- 
vide more data with which to refine the chemical model 
predictions. In particular, the observation of molecules 
such as II2CO and SO, or of the isotope II^'^CO+, would 
help to clarify the nature of the IICO+ ( J=1^0) emission. 
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